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ABSTRACT 


Model  too  to  havo  boon  porformod  to  determine  tho  blast 
v/avo  propagation  in  tho  tunnol  system  of  underground 
□ingle  chnmbor  otorago  sites.  A total  of  nine  configura- 
tions woro  tooted  with  variable  tunnol  diameters  (6,7  - 
13,5  cn),  chambor  volumes  (7250  - 15200  <yo^),  and  TUT 
loading  denoitieo  (1-70 

Analysis  of  the  blast  wave  data  produced  oimplo  scaling 
rolntionohipo  which  contain  a wido  rango  of  important 
goomotrical  paramotora.  Conparioono  arc  givon  with  one 
largo  ocalo  teot  and  tho  importance  of  wall  roughness 
attenuation  is  discussed, 
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1 . INTRODUCTION 

This  report  is  the  third  in  a series  of  five  /I  - 5/ 
describing  the  results  from  an  extensive  series  of 
model  tests  on  underground  ammunition  storage.  In  the 
preceeding  report*  (Report  II  A),  measurements  were  pre- 
sented on  chamber  pressure.  The  present  report  deals 
with  the  problem  of  blast  wave  propagation  from  the 
storage  chamber  down  the  main  passage-way  fcr  typical 
single  chamber  storage  sites. 

The  use  of  the  present  model  tests  in  the  design  of 
underground  ammunition  storage  sites  is  made  possible 
by  two  special  circumstances:  First,  the  models  used 

in  the  tests  we  re  designed  to  represent  a wide  range 
of  typical  full  scale  sites.  Second,  the  models  were 
carefully  constructed  to  have  sufficient  strength  for 
relatively  high  loading  densities  to  cover  a wide  range 
of  full  scale  cases. 

The  fundamental  scaling  relationships  used  in  the  tests 
are  reviewed  briefly  in  Sec.  2,  followed  by  a review 
of  the  experimental  details  in  Sec.  3.  The  principal 
experimental  results  are  presented  in  Sec.  4,  followed 
by  a summary  and  conclusions  in  Sec.  3. 

2.  THEORETICAL  CONSIDERATIONS 

2. 1 Scaling  Relationships 

The  basis  for  the  scaling  relationships  proposed  in  the 
present  report  is  the  familiar  Hopkinson  scaling  laws 
/6/  which  were  reviewed  in  Report  I / 1 / . However, 
for  the  present  experiments  various  empirical  cross 
checks  have  provided  more  general  relationships.  The 
justification  for  these  extensions  will  be  discussed 
next. 

Earlier  tests  in  straight  tubes  and  tunnels  /7/  have 
shown  that  front  pressure  data  (disregarding  wall  rough- 
ness) scale  approximately  as 

p = f(Q/Vt) 


(2.1a) 
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Referring  to  Fig.  2.1,  a),  Q is  the  explosive 
charge  weight  and  the  total  volume  traversed 
by  the  shock  wave  to  the  observation  point. 

Empirical  fits  for  pressures  ranging  from  0,1  - 
100  bar  have  shown  that  the  functional  relation- 
ship of  Eq.(2.1a)  can  be  approximated  by  (smooth 
wall) : 


p a-  24(Q/Vt)0’66 

3 

with  the  units  for  Q and  in  kg  TNT  and  m , 
respectively. 

likewise,  the  date  for  positive  impulse  I and 
positive  duration  t in  straight  tubes  and 
tunnels  have  been  shown  to  scale  approximately 
as  (disregarding  wall  roughness): 

IA/Q  = g(Q/Vt) 

and 

t+A/Q  = h(Q/Vt)  . 

Here,  A is  the  tube  or  tunnel  cross-sectional 
area. 

Detonation  of  explosive  charges  at  the  dead  end 
of  a straight  tube  or  tunnel,  Fig.  2.1,  b,  will 
produce  the  same  scaling  relationships  as  in 
Eqs.  (2.1a)  - (2.1c)  111.  This  geometry  is  of 
course  not  very  different  from  a single  chamber 
storage  site  shown  in  Fig.  2.1,  c. 

The  natural  extension  of  the  scaling  relation- 
ships in  Eqs.  (2.1a)  - (2.1c)  would  therefore  be 


p = f(Q/Vt,  A./Ac) 

(2. Id) 

IA./Q  = g(Q/Vt,  A./Ac) 

(2. 1e) 

t+A./Q  = h(Q/Vt,  A./Ac)  . 

(2. If) 

(2.1b) 

(2.1c) 
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Here,  A./A^  is  the  ratio  between  the  cross- 

t)  ^ 

sectional  areas  of  the  passageway  and  the  chamber, 

and  is  introduced  as  a dimensionless  adjustable 

parameter.  The  ratio  is  partly  justified  on  the 

grounds  that  Eqs. (2. Id)  - (2. If)  correctly  reduce 

to  Eqs.  (2.1a)  - (2.1c)  for  A. /A  = 1.  It  is  not 

J c 

possible  to  proceed  any  further  at  this  point  with- 
out empirical  cross-checks  using  the  actual  experi- 
mental results.  In  Sec.  4 various  attempts  will  be 
made  to  test  specific  functional  relationships  of 
Eqs.  (2. Id)  - (2. If). 

Clearly,  there  are  various  aspects  of  the  blast  wave 
formation  and  propagation  down  the  tunnel  which  are 
not  accounted  for  by  the  simple  scaling  relationships 
proposed  here.  Some  of  these  effects  will  be  dis- 
cussed in  the  following  sections. 


2.2 


Peak  Pressure  Close  to  Storage  Chamber 
Porzel  /8/  has  argued  that  for  a surface  wave  entering 
a tunnel  as  in  Pig.  2.2a,  the  material  velocity  in 
the  direction  of  the  tunnel  is  stagnated  near  the 
entrance.  The  diffraction  into  the  tunnel  is  then 
equivalent  to  the  inverse  of  normal  reflection,  i.e. 
referring  to  Fig.  2.2a,  the  pressure  in  the  tunnel, 

PD,  is  related  to  the  surface  pressure  pc,  through 
the  usual  reflection  formula  (ideal  gas,  y = 1,4): 


2pD(4pD  + 7) 
7 + PD 


(2.2a) 


To  calculate  the  tunnel  diffraction  pressure 
one  needs  the  inverse  relation  to  Eq.  (2.2a): 


PD  " 


14  - 


224p.  - (14  - 


V' 


/ 1 6 (2.2b) 


The  situation  discussed  here  is  similar  to  the 
present  case  where  the  chamber  pressure  is  diffracted 
into  the  tunnel.  This  hypothesis  will  be  tested  in 
Sec. 4. 


A build-up  in  peak  pressure  over  the  diffraction 
pressure  in  Eq.  (2.2b),  is  expected  to  occur  during 
a few  tunnel  diameters  (L/D-  2-3)  just  inside  the 
tunnel  entrance.  This  is  mainly  due  to  multiple  re- 
flections and  the  formation  of  a Mach  stem  from  numerous 
diffracted  blast  waves  /8/.  This  build-up  pressure,  p,, 

Jj 

say,  has  been  determined  empirically  /8/  to  be  about 
40 io  of  the  diffracted  pressure,  i.e. 

PB-  1,4  PD  (2.2c) 

It  should  be  pointed  out  that  the  pressures  discussed 
here  are  average  pressures  corresponding  to  the 
interpretation  of  average  peak  pressures  from  actual 
pressure-time  recordings  / 1 / . 

Beyond  this  build-up  region  in  the  tunnel,  there  is  a 
formation  of  a "turbulent  choke".  Due  to  the  wall 
roughness,  turbulent  surface  layers  are  formed 
which  grow  and  finally  converge  as  shown  in  Fig.  2.2b. 
Y/hen  these  layers  collide,  the  material  flow  at  the 
centre  is  stagnated.  This  is  denoted  a "turbulent 
choke".  As  the  choke  forms  at  a constant  distance 
behind  the  shock  front,  the  flow  is  then  a steady- 
state  phenomenon. 

2 . 3 Effects  of  Wall  Roughness 

The  scaling  relationships  in  Sec.  2.1  obviously  do 
not  account  for  energy  dissipating  effects  due  to 
the  tunnel  walls  such  as 

a)  thermal  energy  transmission  to  the  walls 

b)  elastic  or  inelastic  deformation  of  the  walls 

c)  viscous  loss  due  to  wall  roughness 

This  has  been  discussed  in  Report  I,  and  it  was 
shown  that  it  is  reasonable  to  assume  that  the  energy 
losses  in  a)  and  b)  will  be  relatively  small  compared 
to  the  total  energy  release.  However,  the  loss  due 
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to  wall  roughness  may  significantly  influence 
the  "blast  propagation,  hut  only  over  long  distances 
and/or  in  the  case  of  large  roughness.  A model 
originally  proposed  hy  Porzel  / 8/  states  that  the 
peak  pressure  attenuation  may  he  expressed  as: 

Y(p)  = constant  - 2 £ (e"/D) (L/D)  (2.3a) 

Here,  Y(p)  is  an  impedance  function  defined  as 

Y(p)  = ^ lnp  - | ln(p  + 7)  - £ (2.3b) 

with  p the  peak  pressure  at  a distance  L/D  along 
the  tunnel.  e/D  is  the  average  relative  wall 
roughness,  and  £ is  an  efficiency  factor  in  the 
range  1/2  ^ £ 51,  For  further  details  reference 
is  made  to  the  original  work  /8/  and  Report  I / 1 / . 

For  the  relatively  smooth-walled  steel  tubes  used 
in  the  present  tests,  the  effects  of  wall  roughness 
are  expected  to  he  quite  small.  This  will  he  dis- 
cussed in  more  detail  in  Sec.  4. 

3.  EXPERIMENTAL  DETAILS 

Details  of  the  models,  instrumentation  and  data  re- 
duction are  to  he  found  in  Report  I /I/,  and  we  shall 
therefore  only  summarize  the  main  features  here. 

The  test  programme  was  designed  to  determine  the  blast 
wave  propagation  in  a straight  transport  tunnel  leading 
from  a single  storage  chamber.  A total  of  nine  combina- 
tions of  chamber  volumes  and  tunnel  cross  sections  were 
used  as  shown  in  Fig.  3. a.  Most  loading  densities  varied 

'Z 

between  3 and  70  kg/nr  and  pressure-time  history  was 
recorded  in  most  cases  at  six  different  distances  from 
the  storage  chamber,  ranging  from  approximately  2 to 
80  tunnel  diameters. 

The  pressed  TNT  charges  used  in  the  tests  were  suspended 
in  the  middle  of  the  detonation  chamber  and  initiated 
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with  electrical  blasting  cap  no.  8.  This  has  an 
equivalent  TNT  weight  of  1,5  - 0,5  g.  The  smallest 
charges  used  in  the  tests  we re  25  g TNT  and  the  un- 
certainty in  the  loading  density  stemming  from  the 
blasting  cap  is  therefore  less  than  2$.  However,  the 
total  uncertainty  in  the  loading  density  is  estimated 
to  be  approximately  5$  hue  to  the  spread  in  the 
weights  of  the  standard  TNT  blocks  used. 

To  measure  the  pressure-time  history  in  the  models, 
standard  measurement  techniques  were  used  and  the 
blast  wave  parameters  were  evaluated  using  special 
computer  programmes  as  discussed  in  Report  I / 1 / . 

The  total  uncertainty  in  the  peak  pressure  determina- 
tion was  estimated  to  be  typically  about  10$,  whereas 
the  uncertainty  in  the  impulse  and  positive  duration 
was  estimated  to  be  about  30$.  Pressure-time  records 
obtained  near  the  detonation  chamber  contained  signi- 
ficant ringing  due  multiple  refelcted  blast  waves  and 
the  peak  pressure  determination  was  probably  no  better 
than  approximately  - 20$. 

4.  EXPERIMENTAL  RESULTS  AND  ANALYSIS 

A total  of  approximately  700  blast  wave  parameter 
values  (peak  pressure,  impulse  and  positive  duration) 
were  obtained  for  the  9 configurations  discussed  in 
Sec.  3.  Only  the  most  significant  results  will  be  pre- 
sented here  since  part  B of  this  report  contains  all 
the  pressure-time  recordings  and  tabulated  blast  wave 
parameters . 

4 . 1 Peak  Pressure  Close  to  Storage  Chamber 

Pressure  measurements  were  made  near  the  exit  of  the 
chamber  or  more  specifically  for  L/D  = 2,5  - 2,8  for 
all  the  configurations  listed  in  Pig.  3. a.  These 
results  are  shown  in  Pig.  4.1a  - 4.1c. 

The  dashed  lines  pc  in  these  figures  are  the  average 
peak  chamber  pressures  versus  loading  density 
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which  have  been  presented  in  Report  II  A /2 /. 

The  solid  lines  Pg  are  the  theoretical  prediction  of 
the  diffracted  pressure  using  Eq.(2.2b)  and  the  dotted 
lines  Pg  the  diffracted  build-up  pressure  using  both 
Eqs.(2.2b)  and  (2.2c).  As  may  be  seen  there  is  fair 
agreement  between  the  measurements  and  the  Pg-curves 
considering  the  rough  approximations  made  in  the 
theoretical  predictions  and  the  uncertainties  in- 
volved in  the  interpretation  of  the  average  peak 
pressures  from  measurements.  It  may  be  noted  that 
the  experimental  data  appear  to  be  increasingly 
higher  in  relation  to  the  Pg-curves  for  decreasing 
values  of  A^/Ac.  This  may  be  due  to  an  increased 
build-up  pressure  for  smaller  A^/Ac~values . 

The  same  tendency  appears  to  be  the  case  also  for 
increasing  chamber  volumes.  However,  within  approxi- 
mately - 30 io  the  experimental  data  are  reproduced 
by  the  theoretical  curves. 

4 . 2 Peak  Pressure  Scaling  Relationships 

The  basic  quantity  of  interest  in  the  present 
model  tests  is  the  peak  side-on  overpressure  p in 
the  passageway.  In  Pig.  4.2a,  p is  plotted  versus 
L/D  with  the  loading  density  Q/V^  as  parameter  for 
one  typical  configuration.  Pig.  4.2b  shows  the 
same  data  plotted  versus  Q/Vp  with  L/D  as  para- 
meter. These  plots  are  in  socalled  dimensionless 
form,  i.e.  independent  the  linear  scaling  factor  n. 

(See  Report  I / 1 / ) . 

Letting  subscripts  P and  I!  denote  full  scale 
and  model  respectively, 

Pp  = %,  (4.2a) 

and  Lg/Dg  = nL^/nD^  = I^/D^  • 

For  the  loading  density, 

Qp/^p  = n Q^[/n  = Qjvj/Vjyj 


(4.2b) 
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The  results  in  Figs.  4.2a  and  4.2b  may  there- 
fore in  principle  be  used  for  any  geometrically 
similar  systems.  However,  it  is  of  interest  to 
reduce  the  number  of  curves  necessary  to  characterize 
the  blast  wave  propagation  in  this  particular  system. 
In  Sec.  2.2,  it  was  proposed  to  scale  the  data  accord- 
ing to  Q/Vt,  where  now  is  the  total  volume  of  the 
chamber  and  the  tunnel  out  to  the  observation  point. 
Fig.  4.2c  shows  the  results  in  this  presentation 
for  the  same  data  as  before.  As  may  be  seen,  the 
data  points  are  evenly  distributed  around  one  straight 
line.  This  means  that  one  universal  curve  appears  to 
reproduce  the  tunnel  pressure  at  any  point  in  the 
tunnel  for  any  charge  Q for  the  present  range  of 
measurements.  In  particular,  it  is  remarkable  to 
note  that  even  the  data  for  L/D  = 2,8  in  Fig.  4.2c, 
appear  to  fall  together  with  the  bulk  of  the  data 
for  L/D  > 10. 

Figs.  4. 2d  - 4.2f  shov;  the  variation  of  p versus 
Q/V^.  for  all  the  data  grouped  together  for  the  same 
value  of  A. /A  with  V.  as  parameter.  As  may  be 

J C 1 

seen, there  appear  to  be  no  distinct  systematic 
differences  between  the  three  sets  of  data  for  the 
three  different  chamber  volumes.  However,  there 
is  a systematic  lowering  of  the  bulk  of  the  data 
points  for  decreasing  values  of  A^/Aq.  These 
observations  will  be  tested  more  systematically 
in  the  following  sections  which  contains  the  re- 
sults from  various  empirical  fits  of  all  the  data 
to  the  general  relationship  in  Eq.  (2. Id).  Special 
computer  programs  were  used  for  this  purpose,  and 
for  details  reference  is  made  to  Report  I / 1 / . 

4.2.1  Scaling_Relationshigs_for_Each  Configuration 

This  section  lists  three  methods  examined  for 
scaling  the  peak  pressure  data  for  each  of  the 
nine  configurations  of  the  tests. 
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Method  a 

The  first  method  of  scaling  that  was  investigated 
is  based  on  the  empirical  observation  made  in 
Sec.  4.2.  Here,  it  was  shown  that  for  one  con- 
figuration, log  p varied  approximately  linearily 
with  log  (Q/V^). 

Matematically , this  can  be  expressed  as 

B? 

P = B.,(Q/Vt)  * (4.2.1a) 


where : 

p = measured  peak  pressure 
Q/V^  = effective  loading  density 

= scaling  factor 

B 2 = scaling  exponent 


The  computer  program  made  least  squares  fits  of 
all  the  data  to  Eq.(4.2.1a)  for  each  of  the  nine 

configurations  (A1 , ) , (A1 , V^)f in 

Fig.  3. a to  determine  and 


The  results  of  this  procedure  are  shown  in 
Table  4.2.1a.  As  may  be  seen,  there  appears  to 
be  no  systematic  variations  in  for  different 
chamber  volumes  for  the  same  value  of  A^/Ac. 
However,  on  the  average  there  is  a relatively 
small  increase  in  the  values  of  B£  for  increasing 
chamber  volumes.  For  decreasing  values  of  A. /A  , 
there  is  on  the  average  a distinct  systematic  de- 
crease in  B^  whereas  B ^ increases  somewhat. 


Method  b 

The  second  method  used  to  scale  the  peak  pressures 
is  also  based  on  Eq.  (4.2.1a),  but  including  also 
the  possibility  of  wall  roughness  attenuation 
employing  the  impedance  function  Y(p)  discussed  in 
Sec.  2.3: 


11 


c 


2 


P0  = VQ/Vt) 


(4.2.1b) 


Y(P0)  - Y(p)  = (C3/D)(L/D) 


(4.2.1c) 


where : 


p = fictitions  peak  pressure  without 
wall  friction 

p = measured  peak  pressure • 

Y(p  ) = impedance  function  of  pQ,  Eq.(2.3b) 

Y(p)  = impedance  function  of  p,  Eq.  (2.3b) 

= scaling  factor 

C2  = scaling  exponent 

C 3 = 2 Be  = effective  wall  roughness, 

^ see  Report  I. 


As  in  method  a),  the  data  for  each  of  the  nine  con- 
figurations were  fitted  to  Eqs.(4.2.1b)  and  (4.2.1c). 

The  results  of  this  procedure  are  shown  in  Table  4.2.1b. 
As  may  be  seen,  the  scaling  parameters  and  do 
not  differ  significantly  from  the  earlier  values  for 
and  B^  in  Table  4.2.1a.  The  effective  wall  rough- 
ness is  on  the  average  quite  small  which  is  consist- 
ent with  the  use  of  the  relatively  smooth-walled  steel 
tubes  in  the  tests. 

Method  c 

The  third  method  used  was  identical  to  method  b) 
except  for  introducing  an  additional  parameter  D^: 


The  use  of  is  a crude  way  to  check  the 
possible  effects  of  the  oxygen  deficiency  in 
TNT  / 2/ . 

Again,  the  data  for  each  of  the  nine  configura- 
tions were  fitted  to  Eqs.(4.2.1d)  and  (4.2. 1e). 
The  results  of  this  procedure  are  shown  in 


p0  = VQ/Vt  - V 


(4.2. Id) 


Y(p0)  - Y(p)  = (B3/D)(L/B) 


(4.2. 1e) 
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Table  4.2.1c.  As  may  be  seen,  there  are  significant 
variations  in  the  fitted  values  for  D^,  introducing 
large  variations  in  the  other  scaling  parameters. 

4.2.2  Scaling  Relationships  for  Same  Values  of  A. /A 

J c 

This  section  lists  three  methods  examined  for 
scaling  the  combined  peak  pressure  data  for  V^ , 
and  for  each  value  of  A^/Ac#  This  was  partly 
justified  from  the  findings  in  Sec.  4.2,  which  showed 
that  the  fitted  scaling  parameters  appeared  to  be 
relatively  insensitive  to  the  variation  in  the 
chamber  volume.  The  form  of  the  scaling  relation- 
ships were  identical  to  those  tested  in  Sec.  4.2.1. 

Method  a 

The  use  of  Eq.  (4.2.1a)  produced  the  set  of  scaling 
parameters  shown  in  Table  4.2.2a.  (Note  that  the 
scaling  parameters  now  are  denoted  E^  and  E^ 
instead  of  and  B2).  The  systematic  variation 
of  E.j  and  E2  with  A^/Ac  is  evident. 

Method  b 

The  use  of  Eqs.  (4.2.1b)  and  (4.2.1c)  produced  the 
set  of  scaling  parameters  shown  in  Table  4.2.2b 
(C's  replaced  by  E's).  As  before,  the  effective 
wall  roughness  is  small  as  it  should. 

Method  c 

The  use  of  Eqs.  (4. 2. Id)  and  (4.2. 1e)  produced  the 
set  of  scaling  parameters  shown  in  Table  4.2.2c. 

(D’s  replaced  by  G's).  The  large  variations  in  G^ 
produce  large  scatter  in  the  other  scaling  parameters. 

4.2.3  Scaling  Relationship  for  All  Configurations 

This  section  lists  three  methods  examined  for  scaling 
all  the  peak  pressure  data  at  once  including  also 
the  effects  from  the  variation  of  A ../A  in  the  fits. 
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Method  a 

The  first  method  is  essentially  the  same  as  method 
a)  in  Secs.  4.2.1  and  4.2.2  except  for  introducing 


•The  use  of  is  a way  to  find  an  empirical  de- 


fitting procedure  are  shown  in  Table  4.2.3a.  As 


uncertain. 

Method  b 

The  second  method  is  essentially  the  same  as 
method  b)  in  Secs.  4.2.1  and  4.2.2  except  for 


The  results  of  this  fitting  procedure  are 
shown  in  Table  4.2.3b.  There  is  a consider- 
able improvement  in  the  fitted  parameters  com- 
pared to  method  a). 

Metod  c 

The  third  and  last  method  is  again  essentially 
the  same  as  method  c)  in  Secs.  4.2.1  and  4.2.2 
except  for  introducing  the  A^/Ac-dependence : 


an  additional  parameter  H^: 


H 


5 


(4.2.3a) 


pendence  of  p on  A./A  . The  results  from  this 

J C 


may  be  seen,  the  fitted  value  for  is  quite 


J 


5 


(4.2.3b) 


Y(p0)  - Y(p)  = (J3/D)(L/D) 


(4.2.3c) 


J 


5 


(4.2.3d) 


Y(p0)  - Y(p)  = (J3/D)(L/D) 


(4. 2. 3e) 
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The  result  of  this  fitting  procedure  are  shown 
in  Table  4.2.3c.  The  introduction  of  the  para- 
meter causes  minor  changes  in  the  other  para- 
meters compared  to  method  b) . 

4.2.4  Evaluation  of  the  Empirical  Fits 

On  the  basis  of  the  results  from  the  nine  differ- 
ent empirical  fits  found  in  the  preceeding  sections, 
there  are  three  major  conclusions  which  may  be 
reached: 

i)  The  nine  different  scaling  relationships  all 
reproduce  the  data  within  approximately  twice 
the  expected  experimental  scatter  (-  20 7&). 

There  are  only  minor  differences  for  the  differ- 
ent methods  in  the  data  scatter  around  the 
fitted  curves. 

ii)  Methods  b)  and  c)  involving  the  use  of  the 
proposed  impedance  function,  have  produced 
relatively  small  empirical  values  for  the 
average  wall  roughness.  This  is  consistent 
with  the  relatively  smooth-walled  steel  tubes 
used  in  the  tests. 

iii)  The  introduction  of  an  additional  parameter 
(D4,  G I4)  in  method  c)  to  see  the  possible 
effects  of  the  oxygen  deficiency  in  TNT,  did  not 
improve  the  fits. 

As  a general  conclusion  it  is  therefore  recommended 
to  use  the  results  in  Table  4.2.3b  in  the  prediction 
of  the  smooth  wall  peak  pressures,  i.e.: 

p = (12, 1±0,6)(Q/V. )0,607i°'0l6(A./A  )0» 19-0,04 
o u D c 


• # • 


(4.2.4) 
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Here,  a 2 standard  deviation  (96%  confidence 
interval)  has  been  used  on  the  scaling  para- 
meters. Eq.  (4.2.4)  has  reproduced  all  the 
data  with  an  average  error  of  approximately 
14%.  This  is  quite  satisfactory  considering 
the  estimated  10%  experimental  uncertainty  in 
the  peak  pressure  determination. 

There  are  four  major  limitations  imposed  on  the  present 
results  for  which  the  simole  relationshio  in  Eq . ( ^ . 2 . *0 
is  exnected  to  be  invalid;  a)  different  geometry  of  the 
chamber  exit;  b)  long  chambers;  c)  short  distances  from 
the  chamber  exit  (L/D-*-0);  and  d)  significant  w all 
roughness.  In  regard  to  case  a),  earlier  shock  tube 
experiments  have  shown  that  abrupt  cross-section  changes 
(i.e.  sharp  corners  as  in  the  present  tests)  attenuate 
the  blast  wave  much  more  than  smooth  transitions  79/ . As 
for  case  b),  long  chambers  would  allow  the  blast  wave  to 
travel  in  both  directions  if  the  charge  was  placed  in  the 
middle  of  the  chamber  and  a double-oeaked  blast  wave  would 
be  observed  in  the  tunnel.  One  peak  would  be  due  to  the 
direct  blast  wave  and  the  second  due  to  the  reflected 
wave  off  of  the  back  wall  of  the  chamber.  In  case  c)  the 
complicated  reflection  pattern  of  the  blast  wave  going 
from  the  chamber  to  the  tunnel  precludes  any  unambiguous 
definition  of  a peak  oressure. 

For  significant  tunnel  wall  roughness,  case  d) , 

Eq.  (4.2.4)  may  be  supplemented  by  the  use  of  the 
impedance  function  as  outlined  earlier  in  this 
section  and  in  Report  I / 1 / . As  has  been  discussed 
elsewhere  / 10/,  the  wall  roughness  attenuation  may 
be  considerable  for  typical  tunnels  blasted  out  of 
rock,  and  the  use  of  Eq.  (4.2.4)  will  be  quite  con- 
cervative  or  on  the  safe  side. 
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4. 3 Scaling  Relationships  fcr  Impulse 

The  results  for  the  impulse  are  much  more  difficult 
to  combine  empirically  in  a straightforward  way 
than  those  for  the  front  pressure.  This  is  partly 
due  to  the  significantly  larger  uncertainties  in- 
volved here  and  special  features  are  not  easily 
discernible.  This  is  clear  from  the  .typical  re- 
sults in  Fig.  4.3a.  On  the  average,  the  impulse 
decreases  slowly  with  distance  from  the  storage 
chamber,  but  there  is  no  systematically  large  drop 
near  the  chamber  as  was  observed  for  front  pressure. 

To  test  the  proposed  scaling  hypothesis  in  Sec.  2.1, 

Eq.  (2.1.e),  the  results  for  the  impulse  measure- 
ments are  shown  in  Figs.  4.3b  - 4.30  sis  I A j /Q 
versus  Q/V^.  As  expected,  the  data  show  relatively 
large  scatter,  but  to  a first  approximation  the  re- 
sults fall  on  discrete  curves  with  L/D  as  a para- 
meter. Part,  or  all  of  the  L/D  dependence  may  be 
due  to  the  wall  roughness,  but  so  far  it  has  not  been 
possible  to  account  for  this  theoretically.  The  re- 
sults also  show  that  the  scaled  impulse  is  rela- 
tively insensitive  to  the  values  for  A../A  and  the 
chamber  volume  . As  for  front  pressure,  various 
empirical  scaling  relationships  with  Eq.  (2.1e)  as 
basis  were  fitted  to  the  experimental  results.  In 
particular,  a reasonable  fit  to  all  the  data  v/as 
obtained  using  the  following  expression: 

k2  k3 

I A./Q  = K.,(Q/Vt)  (A,/Ac)  exp(-K4.L/D)  ...(4.3) 

Here,  , P'2  K,,  and  were  introduced  as  scaling  para- 
meters. The  result  of  this  nrocedure  are  shown  in  Table 
4.3.  As  expected,  the  scaling  narameters  are  relatively 
uncertain,  but  Eq.  (4.3)  reproduces  all  the  data  with  an 
average  scatter  of  aonroximatelv  25  %•  This  is  quite 
satisfactory  considering  the  relatively  large  uncertainties 
involved  in  the  interpretation  of  the  impulse. 
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The  simple  relationship  in  Eq.  (4.3)  is  expected  to  be 
invalid  for  the  same  cases  as  listed  in  Sec.  4.2.4  for 
the  peak  pressure  data.  For  significant  wall  roughness  as 
for  tunnels  blasted  out  of  rock,  the  use  of  Eq . (4.3) 
will  probably  produce  result  v/hich  are  quite  conservative 
• or  on  the  safe  side. 


4 . 4 . Scaling  Relati on s h ins  for  Po sitive  Duration  . 

The  results  and  the  interpretation  of  the  positive 
duration  data  show  similar  characteristics  as  for  the 
impulse  data  discussed  in  Sec.  4.3.  The  comments  made 
there  essentially  also  apply  here.  Some  typical  results 


are  shov/n  in  Fig.  4.4.  Again,  there  is  large  scatter 
in  the  data,  but  to  a first  approximation  the  results 


fall  in  discrete  curves  with  L/D  as  parameter.  As  for 
the  impulse  data,  a reasonable  fit  to  all  the  data  was 
obtained  using  an  expression: 


t A./Q  = L 
J 


L L 

(Q/Vt)  * (A./Ac)  exp  (-1^  L/D) 


(4.4) 


Here,  , L^  and  L^  act  as  scaling  parameters.  The 

results  of  the  nonlinear  least  squares  fit  of  Eq.  (4.4) 
to  all  the  data  are  shown  in  Table  4.4.  This  fit  repro- 
duces all  the  results  with  an  average  scatter  of 
approximately  30  % which  again  is  quite  satisfactory 
considering  the  large  uncertainties  involved  in  the 
interpretation  of  the  experimental  results.  The  limita- 
tions are  essentially  the  same  as  those  listed  in  Sec. 
4.2.4  for  the  peak  pressure  data. 


4.5-  Comear ison  with  Large  Scale  Tests . 

Parts  of  some  large  scale  tests  performed  by  the  Nor- 
wegian Defence  Research  Establishment  (NDRE)  /ll/, 
correspond  approximately  to  the  model  tests  reported 

1 1 
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here.  TNT  charges  ranging  from  100  to  1000  kg  were 

positioned  at  a blind  end  of  a rock  tunnel  with  a cross 

2 

section  of  5,5  m . The  pressure  was  recorded  at  6 differen 

distances  from  the  charge  ranging  from  20  to  80  n as 

shown  in  Fig.  *t.5a.  This  situation  would  correspond  to 

the  special  case  where  A. /A  = 1. 

J c 

However,  as  indicated  in  Fig.  iJ^a,  a side  chamber  is 
connected  to  the  tunnel,  and  this  will  ohviously  attenuate 
the  blast  ’wave  somewhat,  but  probably  less  than  20  % / 9 / . 
Fig.  4.5b  shows  peak  pressure  versus  effective  loading 
density  with  the  distance  from  the  charge  in  tunnel 
diameters,  L/D,  as  parameter.  As  may  be  seen,  these 
results  are  consistently  lower  than  the  model  data  for 
A./A^  = 0,^5.  Moreover,  the  large  scale  data  show  a 
pronounced  L/D-denendence  'which  is  undoubtably  due  to  the 
wall  roughness  attenuation  which  was  discussed  earlier. 

Employing  the  impedance  function  discussed  earlier,  the 
large  scale  data  in  Fig.  4.5h  were  fitted  to  Eqs.  (4. 2.1h) 
and  (4.2.1c),  which  produced  a "smooth  wall"  relation- 
ship: 

p0  = (12  - 2)(Q/Vt)0’57  1 °’°4  (4.5a) 

and  an  effective  wall  roughness 

2 £e  = 0, 19  - 0,02  m (4.5b) 

The  corrected  "smooth  wall"  data  p are  shown  in 

o 

Fig.  4.5c  together  v/ith  the  fitted  curve  expressed 
by  Eq.  (4.5a). 

Even  the  "smooth  wall"  data  are  consistently  lower 
than  the  corresponding  model  data,  which  may  be 
due  to  additional  losses  in  the  full  scale  case  • • 

as  large  amounts  of  water  was  present.  This  will 
be  discussed  in  more  detail  in  Report  V in 
connection  v/ith  other  large  scale  tests  in  the 
same  installation. 
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SUMMARY  AND  CONCLUSIONS 

The  methods  in  this  report  for  the  prediction  of 
blast  wave  propagation  in  single  chamber  storage 
sites  constitute  a major  source  of  information  in 
the  planning  of  full  scale  installations  in  rock. 

The  results  provide  a simple  and  sufficiently 
accurate  method  for  the  determination  of  the  blast 
wave  parameters  at  any  point  in  the  transport  tunnel 
for  a wide  range  of  loading  densities  and  geometries. 

This  is  made  possible  by  the  use  of  simple  scaling 
relationships  with  the  following  characteristics: 

a)  All  the  formulas  are  given  in  dimensionless  form 
independent  of  scale. 

b)  The  net  explosive  quantity  (TNT)  is  used  as  basic 
input  parameter. 

c)  The  effects  from  a wide  range  of  geometrical  con- 
figuration (chamber  volumes,  tunnel  cross-sections) 
have  been  introduced  in  a relatively  simple  form. 

d)  The  results  are  valid  for  essentially  smooth-walled 
tunnels  and  are  expected  to  be  conservative  or  on  the 
safe  side  for  tunnels  blasted  out  of  rock. 

e)  The  scaling  relationship  for  the  peak  pressure  data 
may  serve  as  a starting  point  for  the  prediction  of 
the  blast  reducing  effects  of  wall  roughness. 

Direct  comparisons  of  the  present  results  have  been 
made  with  one  large  scale  test  with  a similar  tunnel 
system,  but  with  a ruch  larger  relative  wall  roughness. 
The  relatively  large  difference  in  the  peak  pressure 
attenuation  for  the  model  and  large  scale  tests,  has 
been  partly  accounted  for  as  being  due  to  the  differ- 
ence in  wall  roughness. 

Finally,  it  is  beleived  that  this  study,  with  its  method 
for  the  prediction  of  blast  at  the  tunnel  exit,  is  a 
promising  starting  point  in  the  determination  of  the" 
hazardous  area  around  an  underground  ammunition  storage 
site. 
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Table  4.2.1a 


Least  squares  fits  of  Eq.  (4.2.1a)  to 
the  peak  pressure  data  for  each  of  the 
nine  test  configurations. 


V. 

1 

(cm5) 

VAc 

■D 

B1 

1 

B2 

ErrorB 

(fO  . 

7250 

0,45 

14,4  ± 1,0  | 

0,45  i 0,02 

6 

10900 

0,45 

11,6  ± 1,1  | 

0,51  1 0,02 

11 

15200 

0,45 

10,5  - 1,1  i 

0,62  - 0,02 

12 

7250 

0,25 

10,0  - 1,1 

0,52  - 0,02 

9 

10900 

0,25 

7,7  i 1,1 

0,64  - 0,02 

10 

15200 

0,25 

8,5  - 1,1 

0,66  - 0,02 

10 

7250 

0,11 

6,9  - 1,1 

0,64  - 0,02 

10 

10900 

0,11 

6,7  - 1,1 

0,71  - 0,02 

8 

15200 

0,11 

7,1  i 1,1 

0,70  - 0,02 

13 

In  units  where  p is  expressed  in  bar,  Q in  kg  TNT, 

3 

and  in  m . 

^ Average  difference  between  fitted  curve  and  experi- 
mental data. 


Table  4.2.1b. 


Nonlinear  least  squares  fits  of  Eqs.  (4.2.1b) 
and  (4.2.1c)  to  the  peak  pressure  data  for 
each  of  the  nine  test  configurations. 


vi  AA 

(cm5) 

a 

c. 

I 

1 

1 

1 

1 

1 

1 

1 

1 C\J 

1 o 

I 

1 

l 

1 

1 

C3 

(mm) 

T) 

Error 

(*) 

7250  0,45 

11,7  - 1,2 

0,50  - 0,04 

•s 

o 

+ 1 

o 

•» 

o 

12 

10900  0,45 

10,0  - 0,9 

0,58  ± 0,03 

0,1  i 0,3 

15 

15200  0,45 

10,0  t 0,7 

0,63  - 0,02 

0,0  - 0,3 

13 

7250  0,25 

10,7  1 0,7 

0,51  - 0,02 

0,22  - 0,14 

9 

10900  0,23 

8,4  - 0,7 

0,63  - 0,02 

0,17  - 0,12 

9 

15200  0,23 

9,6  ± 0,6 

0,64  - 0,02 

0,31  1 0,11 

9 

7250  0,11 

6,9  - 0,6 

0,63  - 0,02 

0,00  - 0,07 

10 

10900  0,11 

7,0  - 0,8 

0,71  1 0,03 

0,11  - 0,09 

8 

15200  0,11 

7,8  ± 0,6 

0,680i  0,018 

0,13  1 0,08 

13 

a '5 

In  units  where  p is  expressed  in  bar,  Q in  kg,  and  V+  in  m . 


b 


Average  difference  between  fitted  curve  and  experimental 
data. 


Table  4.2.1c.  Nonlinear  least  squares  fits  of  Eqs.(4.2.1d)  and 

(4.2. 1e)  to  the  peak  pressure  data  for  each  of 
the  nine  test  configurations. 


a In  units  where  p is  expressed  in  bar,  Q in  kg  TNT,  and  V, 

3 15 

m m . 


b 


Average  difference  between  fitted  curve  and  experimental  data. 


Table  4.2.2a.  least  squares  fits  of  Eq.  (4.2.1a)  to 

the  combined  peak  pressure  data  for  all 
three  chamber  volumes  for  each  value  of 

W 


A. /A 
3 c 

a 

ei 

E2 

" ' b 

Error 

w 

0,45 

11,6  - 1,0 

0,535  1 0,015 

12 

0,25 

8,8  - 1,0 

0,601  - 0,015 

13 

0,11 

6,9  ± 1,0 

0,680  - 0,014 

12 

a 


In  units  where  p is  expressed  in  bar,  Q in  kg  TNT, 
•3 

and  in  m . 


Average  difference  between  fitted  curve  and  experi- 
mental data. 


Table  4.2.2b.  Nonlinear  least  squares  fits  of 

Eas.  (4.2.1b)  and  (4.2.1c)  to  the  com- 
bined peak  pressure  data  for  all  three 

chamber  volumes  for  each  value  of  A./A  . 

J c 


A./A 
3 c 

11 

II 

ll 

11 

Hrj  ll 
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ll 

II 

ll 

II 

It 

It 

II 

11 

II 

11  CvJ 

ll  Ph 

ii 

ii 

it 

ii 

u 

(mm) 

Error 

__(*).  - .... 

0,45 

11,3  1 0,5 

0,545  ± 0,016 

0,00  - 0,16 

13 

0,23 

9,7  - 0,5 

0,588  - 0,015 

0,22  - 0,09 

12 

0,11 

7,3  1 0,4 

0,667  - 0,014 

0,06  - 0,05 

12 

====== 

ii 

II 

II 

II 

11 

II 

11 

11 

11 

II 

11 

ii 

II 

II 

II 

ll 

ll 
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It 
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11 
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II 
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II 

ll 

11 

11 

ll 

ll 

a In  units  where  p is  expressed  in  bar,  Q in  kg  TNT, 

3 

and  in  m . 


b 


Average  difference  between  fitted  curve  and  experi- 
mental data. 


Table  4.2.2c.  Nonlinear  least  squares  fits  of 

Eqs.  (4.2.1a)  and  (4.2.1e)  to  the  com- 
bined peak  pressure  data  for  an<^ 

Vy  for  each  value  of  A^/Ac. 


VAC 

11 

It 

1! 

1! 

1! 

Q II 

it 

ll 

ll 

ji 

a 

G2 

G3 

(mm) 

G4 

Ocg/m3) 

0 

Error 

(?S) 

0,45 

10,3  1 0,9 

0,57  - 0,03 

0,0  - 0,2 

0,14  ± 0,06 

13 

0,23 

9,9  1 0,8 

0,58  - 0,02 

0,22  - 0,09 

0,03  - 0,08 

12 

0,11 

5,1  ± 0,8 

0,78  - 0,03 

0,06  - 0,05 

-0,77  - 0,14 

12 

a 


In  units  where  p is  expressed  in  bar,  Q in  kg  TNT, 
'•5 

and  in  nr. 


b 


Average  difference  between  fitted  curve  and  experi- 
mental data. 


CV, 


Table  4.2.3a.  Nonlinear  least  squares  fits  of 


Eq.  (4.2.3a)  to  the  combined  peak 
pressure  data  for  all  the  configura- 
tions. The  average  difference  between 
the  empirical  curves  and  experimental 
data  is  19/®. 

Parameter 

Numerical  value 

a 

7,6  - 0,8 

H2 

0,68  - 0,03 

It 

II 

It  \J1 

II 

II 

II 

ll 

11 

II 

II 

ll 

tl 

0,03  1 0,04 

a In  units 

where  p is  expressed  in  bar,  Q in 

kg  TNT  and  in  m^. 


a 


Table  4.2.3b. 


Nonlinear  least  squares  fit  of 
Eas.  (4.2.3b)  and  (4.2.3c)  to  the  com- 
bined peak  pressure  data  of  all  the  con- 
figurations. The  average  difference  be- 
tween the  fitted  curve  and  experimental 
data  is  14%. 


Parameter 

Numerical  value 

a 

n 

12,1  i 0,3 

*2 

0,607  ~ 0,008 

I3 

0,00  - 0,05  (mm) 

1 M 
1 VJ1 
I 
! 

1 

1 

I 

I 

1 

l 

1 

u 

0,19  1 0,02 

a In  units  where  p is  expressed  in  bars,  Q in 
kg  TNT,  and  in  nr*. 


Table  4.2.3c. 


Nonlinear  least  squares  fit  of 
Eos.  (4.2.3d)  and  (4.2.3e)  to  the  com- 
bined peak  pressure  data  of  all  the  con- 
figurations. The  average  difference  be- 
tween the  fitted  curves  and  experimental 
data  is  14$. 


Parameter 

Numerical  value 

a 

11,5  1 0,5 

j2 

0,624  * 0,014 

J3 

0,00  - 0,05  (mm) 

h 

-0,08  ± 0,05  (kg/m5) 

J5 

0,19  i 0,02 

In  units  where  p is  in  bar,  Q in  kg  TNT,  and 
•5 

V.  in  nr  . 


a 


Table  4.3. 


Nonlinear  least  squares  fit  of  Eq.  (4.3) 
to  the  combined  impulse  data  of  all  the 
configurations.  The  average  difference 
between  the  fitted  curve  and  experiment- 
al data  is  25^. 


Parameter 

Numerical  value 

a 

ki 

30,2  - 1,5 

k2 

- 0,33  - 0,02 

K3 

- 0,02  - 0,04 

1 M 
i 

t! 

0,0160  - 0,0011 

In  units  where  I is  in  bar  ms,  Q in  kg  TNT, 

2 3 

A.  in  m and  V.  in  m . 

J ^ 


a 


Table  4.4. 


Nonlinear  least  squares  fit  of  Eq.  (4.4) 
to  the  combined  positive  duration  data  o 
all  the  configurations.  The  average  dif 
ence  between  the  fitted  curve  and  experi- 
mental data  is  31$. 


Parameter 

Numerical  value 

a 

h 

20,3  - 1,1 

h 

- 0,81  - 0,03 

b 

0,16  - 0,03 

h 

0,0211  - 0,0009 

In  units  where  t is  in  ms,  Q in  kg  TNT,  A.  in 

2 3 + J 

m and  in  m . 
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Pig.  2.1.  Propoccd  extension  of  coaling  relationships  for 

^ j — V,  . . w'  V,  _ ^ ^ 0 ^ ^ J i.11  Cy  - -•  ^ ^ 

relationships  in  straight  tunnol3,  a)  and  b). 


p , Incident  wave 

^ p Reflected 
C A | wave 

&"  - - ■ 


D 


Fig.  2.2a  Diffraction  of  a surface  blast 
wave  entering  a tunnel  /8/. 
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Fig.  2.2b  Blast  wave  propagation  in  a tunnel 
with  the  formation  of  a turbulent 
choke  /8/. 


V: 


Q 


Ar  -314  cm 


;x^*  -C 

I -Ai  ! 


-t>! 


i 

V. 

l 
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P 

V 

(mb 

0 

A. 

0 

(cm2) 

A/ 

(m2) 

n 

i 

7250 
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D.a. 


Model  configurations  used  in  the  single  chamber 

storage  tests.  The  three  chamber  volumes 

V.  and  the  three  tunnel  cross-sections  A.  given  in 
1 J 'J  " 

the  table  thus  formed  nine  combinations  of  Y.  A..  v:j.th 

J 

i = 1,  2,  3 and  j = 1,  2, 


Choosing  A.  - 20  m ' as  a typical  full  scale  value 

J 

for  the  tunnel,  cross-sec lion,  the  linear  scaling 
factors  n have  the  values  shovm  in  the  last  column . 
The  corresponding  full  scale  chamber  volumes  are 


p (bar) 


Pig.  4.1a  Peak  pressure  near  the  chamber  exit  versus 
loading  density  for  three  chamber  volumes. 

The  ratio  between  the  tunnel  and  chamber  cross- 

sections  is  A. /A  = 0,45. 

J C 

The  dashed  line  represents  the  average  peak  chamber 

pressure,  p , determined  earlier  / 2/ , and  the  solid 
c 

line,  p-r.  is  determined  from  the  inverse  reflection 
formula  in  Eq.(2.2b).  The  dotted  line  is  the 
diffacted  build-up  pressure,  p-o,  using  both 
Eqs.  (2.2b)  and  (2.2c). 


p (bar) 


T 
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Fig.  4.1b  Peak  pressure  near  chamber  exit  versus  loading 
density  for  A. /A  = 0,23  (See  Fig.  4.1a  for  de- 

tails) . 


p (bar) 


Fig.  4.1c  Peak  pressure  near  chamber  exit  versus  loading 

density  for  A. /A  = 0,23  (See  Fig  4.1a  for  de- 
J c 

tails) . 
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i’ig.  4*2a  Peak  pressure  versus  cistunce  from  chaniber 
exit  for  various  loading  densities. 
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dons 
same 


pressure  versus  chamber  loading 
test  configuration  as  in  Fig*  4.2a. 


Fig.  4.2b 
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Fig*  4.2c.  Peal-:  pressure  versus  effective  loading  density 
for  the  3 i.vv^G  test  configuration  as  in  Pig.  4.2a 
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are  also 


/2 / e.  *.d  in  s vc-  A <*;ht  tahas  / 7 / 
uied  i n tiio  figure. 


Pig.  4.2 e 


Peak  pressure  versus  effective  loading  density 

for  A ./A,  = 0,23  (Joe  Fig.  4. 2d  for  details). 

0 e 


p (bar) 


Pig,  4.2f.  Peak  pressure  versus  effective  loading  density 
for  A j/Ac=0, 1 1 . (See  Pig.  4. 2d  for  details). 
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Pi g.  4. 3a. Example  of  positive  impulse  verrtus  distance  fron 
ch umber  exit  expressed  in  turnicO.  aicenetc  rs. 
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Pig.  4.3b  Scaled  impulse  versus  effective  loading  density 
[Eq.(2.1e)3  at  various  distances  in  tunnel  dia- 
meters from  the  chamber  exit. 
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Pig.  4.3c  Sealed  impulse  versus  effective  loading  density 
lEq.(2.1e)J  ax  various  distances  in  tunnel  dia- 
meters from  the  chamber  exit. 
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Fig.  4.3d  Scaled  impulse  versus  effective  loading  density 
Eq.(2.le)  at  various  distances  in  tunnel  dia- 
meters from  the  chamber  exit. 
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Fig.  4.3e  Scaled  imnulse  versus  effective  loading  density 
[,Fq.(2.1e)3  at  various  distances  in  tunnel  dia- 
meters from  the  chamber  exit. 
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Fig.  4.3f  Scaled  iinpulse  versus  effective  loading  density 
[3q.(2.ie)]  at  various  distances  in  tunnel  dia- 
meters from  the  chamber  exit. 
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Fig.  4.3g  Scaled  impulse  versus  effective  loading  density 
[Eq. (2. 1 e)j  at  various  distances  in  tunnel  dia- 
meters from  the  chamber  exit. 
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Fig.  4.3i  Scaled  impulse  versus  effective  loading  density 
[iSq.  (2. 1 e)]  at  various  distances  in  tunnel  dia- 
meters from  the  chamber  exit. 
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Pig.  4. Scaled  impulse  versus  effective  loading  density 
[pq.(2.1e)3  at  various  distances  in  runnel  dia- 
meters from  the  chamber  exit. 
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Fig.  4.5a 


Tunnel  nvuteir.  in  the  ocaio 
formed  by  MDRE  /1 5/. 
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j'ig.  4.5b 


Pea’;  pressure  versus  effective  loading  density 

for  the  large  scale  tests  performed  by  KDR1  > ",  5/ 

in  comparison  v/ith  the  model  data  for  A. /A  - .0,  if. 

J c 

The  L/D-dependence  is  consistent  with  an  average 
wall  roughness  cf  approximately  0,1  r; . 
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p (bar) 
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?i£.  4.5c. 


“Smooth  wall  « peak  pressure  r>  versus  effectiv. 
loading  density  for  the  larce°ocale  test  data 


in  jf'xg . 'r.yc, 
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t 


